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Abstract 

In this paper, we prove that there exists a unique weak solution to the mixed boundary 
value problem for a general class of semilinear second order elliptic partial differential 
equations with singular coefficients. Our approach is probabilistic. The theory of Dirichlet 
forms and backward stochastic differential equations with singular coefficients and infinite 
horizon plays a crucial role. 
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1 Introduction 



In this paper, our aim is to use probabilistic methods to solve the mixed boundary value 
problem for semilinear second order elliptic partial differential equations (called PDEs for 
short) of the following form: 



{Lu{x) = —F{x,u{x),'Vu{x)), on D 
l^{x) - B ■ n{x)u{x) = <^{x) ondD 

The elliptic operator L is given by : 

L = ■{AV) + B-V -V ■{B-) + Q (1.2) 



\ E ^ [a,ix)±-y±B.-Xx)^ d^viB.) + Qix) 
i,i=i * ^ J i=i ' 



on a d-dimensional smooth bounded Euclidean domain D. 

A{x) = {aij)i<ij<d- — >■ R*^ ® is a smooth, symmetric matrix-valued function which is 
uniformly elliptic. That is, there is a constant A > 1 such that 

jidxd < A{-) < Aldxd. (1.3) 

Here B = {Bi,...,Bd) and B = {Bi,...,Bd) : ^ are Borel measurable functions, which 
could be singular, and Q is a real- valued Borel measurable function defined on R^ such that, 
for some p > f , 

Ij,{\B\'' + \B\'' + \Q\)eLP{D). 
L is rigorously determined by the following quadratic form: 

1 v-^ [ du dv -sr^ /" n / \ 

Q{u,v):={-Lu,v)L2i^o) = - aij{x) ——dx - Bi{x)—v{x)dx 

— \^ / Bi{x)^—u{x)dx— I Q{x)u{x)v{x)dx. 
Details about the operator L can be found in [S], [TB] and 



The function F(-, •, •) in (jl.ip is a nonlinear function defined on R^ x R x R'^ and <I>(x) is a 
bounded measurable function defined on the boundary dD and 7 = An, where n denotes the 
inward normal vector field defined on the boundary dD. 

To solve the problem (jl.ip . it turns out that we need to establish the existence and unique- 
ness of solutions of backward stochastic differential equations (BSDEs) with singular coeffi- 
cients and infinite horizon, which is of independent interest. 

Probabilistic approaches to boundary value problem of second order differential operators 
have been adopted by many authors and the earliest work went back as early as 1944 in |12j . 
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There has been a lot of study on the Dirichlet boundary problem (see [T], [8], [3], [6], [TT] and 
|22j). However, there are not many articles on the probabilistic approaches to the Neumann 
boundary problem. 

When A = I, B = and B = 0, the following Neumann boundary problem 

^ Ati(x) + qu{x) =0, on D 
-2^i^)=m ondD 

was solved in [1] and [11], which also gives the solution the following representation: 

roo J 

u{x) = E^[ efo'^(^-)^-^{Bt)dL% 
Jo 

where {Bt)t>o is the reflecting Brownian motion on the domain D associated with the infinites- 
imal generator 

and L^, t > is the boundary local time satisfying = lQ£){Bs)dL'^. 

But when B ^ 0, the term V-(-B-) is just a formal way of writing because the divergence does 
not exist as IB is only a measurable vector field. It should be interpreted in the distributional 
sense. For this reason, the term V • (B-) can not be handled by Girsanov transform or Feyman- 
Kac transform. 

The study of the boundary value problems for the general operator L in the PDE literature 
(see e.g. [9], |20] ) was always carried out under the extra condition: 

-div{B) + Q{x) < 

in the sense of distribution in order to use the maximum principle. 

When F = 0, i.e. the linear case, problem (II. ip was studied in [4]( see also [3j for the 
Dirichlet boundary problem). The term V • {B-) is tackled using the time-reversal of Girsanov 
transform of the symmetric reflecting diffusion {Q.,P^,Xf,t > 0) associated with the operator 

Lo = ^V-{AV). 

The semigroup St associated with the operator L has the following representation (see [5]): 

Stfix) = exp( [\A-'Bnx',)dM^ + ( (\a-^ Bf {X',)dM',) o 7? 

Jo Jo 

-\ [\B-B)A-\B-BnX^)ds+ fQiX^)ds)], 
^ Jo Jo 

where is the martingale part of the diffusion X^ and 7^ is the reverse operator. 
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The main purpose of this paper is to study the nonhnear equation (jl.ip fi.e. F ^ 0), which 
can not be handled by the methods used for the hnear case. Our approach is first to solve a 
backward stochastic differential equation (BSDE) with singular coefficients and infinite horizon 
to produce a candidate for the solution of the boundary value problem and then to show that 
the candidate is indeed a solution. The results we obtained for BSDEs with infinite horizon 
are of independent interest. 

We would like to mention that the first results on BSDEs and probabilistic interpretation 
of solutions of semilinear parabolic PDEs via BSDEs were obtained by Peng and pardoux in 
[19\ . (17] and [18]. There the operator L is smooth and the solution is a viscosity solution. We 
stress that the solutions we considered for PDEs in this paper are Soblev (also called weak) 
solutions, not viscosity solutions. 

In [22] . the corresponding Dirichlet problem for the semilinear elliptic PDEs: 

( Lu{x) = —F{x,u{x),Vu{x)), on D (14) 

was solved. The strategy in [3], [22] is to transform the general operator L by a kind of h- 
transform to an operator of the form: L2 = |V(ylV) + 6 • V + g which does not have the "bad" 
term such as V{B-). This idea is used in current paper too. 

The BSDEs we studied are inspired by the ones in [1^ where the author gave a probabilistic 
interpretation of the solution to the following Neumann problem: 

f (iA - u)u{x) = 0, onD 



The content of the paper as follows. In Section 2, we study the following BSDEs with 
infinite horizon: 

dY{t) = -F{X{t),Y{t), Z{t))dt + e/o '?(^(«))'^*$(X(s))(iLt + {Z {t) , dM {t)) , 

lim e/o'^(^W)'^«y, = in L^fl), (1.5) 

where {X(t))t>o is the reflecting diffusion associated with an infinitesimal generator of the 
form: A = ^V(AV) + 6 • V, M{t) is the martingale part of X{t), Lt is the boundary local 
time of X and d{-) is an appropriate measurable function. The existence and uniqueness of an 
L^-solution (y, Z) is obtained. 
In Section 3, we solve the linear PDEs of the form: 

J lV{AVu){x) + b-\7u{x) + qu{x) = F{x), on D 

[ ^^(x) = <t>{x) ondD. ^^-^f 

under the condition: 

dLt] < 00 

Jo 
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for some xq £ D. Useful estimates for local time and Girsanov density are proved which will 
also be used in subsequent sections. 

In Section 4, we obtain the solution of the semilinear PDE: 

J ^\/{AVu){x) +b-\/u{x) + qu{x) = G{x,u{x),Vu{x)), on D 
\h^i^)=4>{x) ondD. ^^-^f 

To this end, we first use the solution {Yx{t), Zx{t)) of the BSDE (|1.5p to produce a candidate 
uo{x) = Ex[Yx{0)] and then find a solution u of an equation like (II. 6p with a given F{x) := 
G{x, uo{x),vo{x)). Finally we identify u with uq. In Section 5, we consider the general problem: 

Lu{x) = —F{x,u{x)), on D 

^^{x) - B ■ n{x)u{x) = ^{x) ondD ' 

We apply the transformation introduced in [3] to transform the problem (ll.Sp to a problem 
like (jl.7p . An inverse transformation will yield the solution of the problem (jl.Sp under the 
condition that the norm of B is sufficiently small. 

To remove some of the restrictions imposed on B in Section 5, in Section 6, we study the 
L^-solutions of the BSDEs (jl.Sp under appropriate conditions. Our approach is inspired by 
the one in [2]. The study of L^-solutions and L^-solutions of the BSDEs ([LS]) are carried out 
in Section 2 and Section 6 separately because the methods used for these two cases are quite 
different. 



2 BSDEs with Singular Coefficients and Infinity Horizon 

Consider the operator 




on the domian D equipped with the Neumann boundary condition: 

d 

— := (An, V-) = 0, on dD. 

By [13], there exists a unique reflecting diffusion process denoted by (H, J-t, Xx{t), Px,6t, x G D) 
associated with the generator Li. 

Here ^ : — )• is the shift operator defined as follows: 

Xx{s){ef) = Xx{t + s), s,t>o. 

Let Ex denote the expectation under the measure Px- 

Set 6 = bd}, where h = ^ Y.j ^ + ^i- 

Then the process Xx{t) has the following decomposition: 

Xx{t) = Xx{Q) + Mx{t)+ I b{Xx{s))ds+ I An{Xx{s))dL,, Px-a.s.. (2.1) 

JO JO 
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Here Mx{t) is a J-'t square integrable continuous martingale additive functional. And Lt is a 
positive increasing continuous additive functional satisfying Lt = I^X:,{s)edD}dLs- 
We write Xx{t) as X{t) for short in the following discussion. 

In this section, we will study the backward stochastic differential equations with singular 
coefficients and infinite horizon associated with the martingale part Mx{t) and the local time 
Lt- A unique solution of such BSDEs is obtained. 

Let g{uj,t,y,z) : n x R+ x R x R<^ ^ R he a progressively measurable function. Consider 
the following conditions: 

(A.l) {yi - y2){g{t,yi,z) - g{t,y2,z)) < -ai{t)\yi - y2p, 
(A.2) \g{t,y,zi) - g{t,y,Z2)\ < a2\zi - Z2\, 
(A.3) \g{t,y,z)\ < \git,0,0)\ + ^sim + \y\). 

Here ai{t) and 03 (t) are two progressively measurable processes and 02 is a constant. 

Set a{t) = —ai{t) + 60^, for some constant 5 > where A is the constant appeared in (jl.3p . 



Lemma 2.1 Assume the conditions (A.1)-(A.3) and 

roo ^ 

Ex[ e2/o«(«)^«|5(t,0,0)|2dt] 



Then there exists a unique solution {Yx{t), Zx{t)) to the following backward stochastic differ- 
ential equation: 

Yx{t) = Yx{T)+ [ g{s,Yx{s),Zxis))ds- [ < Z,(s), dM,(s) >, t<T; 
Jt Jt 

lim e^*"(")'^"K,(t) = 0, in L^{^). (2.2) 

t— >-oo 

Moreover, 

S^.[supe2i'o'*H'^«|y^(t)|2] < 00 and E^i e^/o '*H'^«|Z^.(s)|2(is] < 00. (2.3) 
t Jo 

Proof. 
Existence: 

The proof of this lemma is similar to that of Theorem 3.2 in [22], but the terminal condi- 
tions here are different. By Theorem 3.1 in [22], the following BSDE has a unique solution 
(Kr(t),Z«(t)): 

g{s,Y^{s),Z2{s))ds- <Zl{s),dMx{s)>, t < n; (2.4) 

and moreover, 

y,"(t) = o, z:(i) = o, t>n. 
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Fix t > and n > m> t. It follows that 

/"n 
Jm 

-2j\^K<-)'^-{Y^[s) - y-(s)) < Z2{s) - Z^{s),dM^{t) > 
Choose two positive numbers Si and 82 such that Si > and di + 82 < S. Then from 

/n 
^2f^ aiu)du^Yi,\s) - Y:\s)){g{s, Y^is), Zl{s)) - g{s, Y^^is), Z:\s)))ds 

/n 
01(5)6^/0 «(")'^"|y;^(s) - Y^{s)\'^ds 

/n 
^2S^a{u)du^Y^^s) - Y^{s)\^ds 

+^J" ^'^°''^''^'''{MX{s)){Z^is) - Z^{s)), (Z^^is) - Z:\s)))ds 

and 

pn 

2 / e2/oa(«)'^«(y«(s)-y-(s))5(s,0,0)ds 

Jm 

< 2S2al /"e^/o «H<i«|yj»(s) - Y^{s)fds + — ^-^ /" e^/o q, 0)pds, 

Jm 2()2a2 Jm 

it follows that 

. -1-1 roo 

E,[e'fo aiu)du\Y-{t) - y-(i)p] + -(1 - I e^/o «H<^«|Z,"(5) - Z:^{s)\'ds] 

^02a2 Jm 

This implies that 

roo 

E^[ e2i'o«(")'^"|Z^(s)-Z^(s)|2ds] ^0, as m,n^ 00. 
Jo 

Hence there exists Z^ such that 

Z^ = \[m e-fo^(^)d^Z^ in L^ffO, 00) x O). 

n— >-oo 



At the same time, we also obtain the following estimates: 

t 

^0202 Jm 

+2sup| / e^/o - Y^{s)) < Z^{s) - Z^{s),dM^{t) > \. 

t Jt 

Taking expectation on both sides of the above inequality, by BDG inequality, we obtain 
t 

1 

/n 

1 /-n It 

2d2a| J„ 2 t 

/•oo 

Jo 



Thus 



E^[supe2~''o - 



-| pn poo 

—E^[ / e^io' «(")<^" 1^(5, 0,0) pds] + 2C2^a;[ / e^/o «H<i«|Z^(s) - Z^{s)\^ds] 

Jm Jo 



— > 0, as m,n ^ oo. 



So, there exists {Yx{t)} such that 

lim E^[sup \%{t) - eio* = 0. 

n— >oo J 

For any e > 0, there exist a positive number such that for any n> N, 



E,[sup fyt) - e/o < £ 



For t> N, noticing Y^{t) = 0, it follows that 

< 2^^[sup \Y^{t) - e/o «M<^«|yjV(^)|2] ^ 2E^[e2/o «W<^«|yjV(^^|2] 

+ 



< e. 
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Thus we have liinE^[\Y^{t)\^] = 0. 



By chain rule, it is easy to see from ()2.4p that 

Y^{t) = e~ ^0 and Z^{t) = e~ 

satisfy the equation ()2.2p and 

hm i?,[e2/o -(")'^''|y,(t)|2] = hm E^[\Y,{t)\'] = 0. 



From the above proof, we also see that (j2.3p holds. 
Uniqueness: 

Suppose that (Y^,Z].) and (y^^,Z^) are two solutions of the equation ()2.2p . 
Set y;(t) = yi(t) - Y^{t) and ^^.(t) = Zl{t) - Zl{t) . Then 

+ ei'o«(")^"(Z,(t),ciM,(t)). (2.5) 

By Ito's formula, we get, for any t < T, 

rT 



-2 a(s)e2/o"(«)'^«|y^(s)|2(is 

-2 j\{s)e''fo-(-)<i-Y,{s){Z,is),dM,{s)) (2.6) 
By condition (A.l) and (A. 2), we have 

2^''e2/o-W'^«y.(.)(5(.,y,i(.),zi(.))-<7(.,y,2(.),z2(s)))ds 

= 2 jJ"^e2/.>H'^"y,(.)(5(.,y,i(.),Zi(.)) -5(.,y,2(.),Zi(.)))(i. 

+ 2 jj'^e^/o '^W'^«y.(.)(5(.,y,2(,),zi(.)) -5(.,y,2(.),z2(.)))d. 

< -2 rai(s)e2/o'^(")'^"|y,(s)|2ds + a2 Te^/o '^H'^"y,(5)|Z,(s)|ds 

Jt Jt 

< -2 j\i{s)e'^fo <^)<i^ 1 1; (s) 1 2rfs + c'a2 W'^" | (s) | ^ds 

+02^ re2/o«W'^«|Z,(s)|2ds. (2.7) 
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Choosing d = 26a2, we obtain 

< e2/o^'^W'^"|y,(r)|2-2 ra(s)e2/o'^H'^-y,(s)(Z,(s),dM,(s) (2. 



Taking expectation on both sides of the above inequaUty, we get that, for any t < T, 
For both and Y^ satisfy the terminal condition in (j2.2p . so that 



r-i>oo 



which leads to ^^[e^ /o "H'i«|y^(t)|2] = o. 

We conclude that Y^{t) = Y^{t) and Z]^{t) = Zl{t). □ 



We now want to apply Lemma l2.1l to a particular situation. 
Let F{x, y, z) : R'^ x R x R'^ R he a Borel measurable function. Consider the following 
conditions: 

(D.l) (yi -y2)iF{x,yi,z) - F{x,y2,z)) < -di{x)\yi -y2?, 
(D.2) \F{x,y,zi) - F{x,y,Z2)\ < ^2^1 - 2^2!, 
(D.3) \F{x,y,z)\ < \F{x,0,z)\+K{x){l + \y\). 
Set d{x) = —di{x) + (5^2 for some constant 5 > 
The follows result follows from Lemma |2.1[ 

Lemma 2.2 Assume the conditions (D.1)-(D.3) and 

/■oo J 

e2/o°'(^("))'^''|F(X(t),0,0)|2dt] <oo. 



Then there exists a unique solution (Y^it), Zx{t)) to the following equation: 

Y,{t)=Y,{T)+ [ F{X{s),Y,is),Z,{s))ds- [ < Z^{s),dM^{s) >, t<T- 
Jt Jt 

lim ei'o'^(^("))'^"y^(t) = 0, in L'^{n). (2.9) 

Consider the following condition instead of {D.3). 
(D.3)' \F{X{t),y,z)\ < K{t), for any y e R and z G R'^. 

Let <I? be a bounded measurable function defined on dD, and function q G LP{D), for p > f • 
The following theorem is the main result in this section. 
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Theorem 2.1 Assume the conditions (D.l), (D.2) and (D.3)\ 

/•CO 

Jo 

for some xq G D and for x £ D, 

/ e^/o dix{u))du^^2j^ q{x{u))du ^ \K{t)\'^}dt] < oo. (2.10) 
Jo 

Then there exists a unique solution (Yx,Zx) to the following BSDE: 

Y,{t) = Y,{T)+ r F{Xis),Y^{s),Z,{s))ds- r e^o^^^(-))'i^^X{s))dL, 
Jt Jt 

{Z,{s),dM,{s)), for t<T, (2.11) 

and 

lim e/o'='(^("))'='"yt = in L^{Q). (2.12) 

Proof. 
Uniqueness: 

Suppose that {Y^,Z^) and (Y^,Z'^) are two solutions of the equation (|2.11|) satisfying (j2.12p . 
Set y^(t) = Y^{t) - Y^{t) and Z^{t) = Zl{t) - Z^{t) . Then 

d{efod(x(^))d'^Y^{t)) = -efod^^(^))d'^{F{X{t),Y^{t),Zl{t))-F{X{t),Y^{t),Z^{t)))dt 
+ d{X{t))e^od{x(u))duY^^^^dt 

+ e^o'i'-^(^))<^^{Z^{t),dM^{t)). (2.13) 
By Ito's formula, we get, for any t < T, 

e'^Io<i(X(^))d^\%(t)\'^ + j\'^Io<X(^))d^{A{X{s))Z^{s),Z^^^^^ 
+ 2 l\'Io^(''(-))^^Y,(s){F{X{s),Y,\s),zU^^^ 

i-T 

- 2 d(X(s))e2/o'^W"))'='"|y^(s)|2ds 

- 2 diXis))e^ /o ^(^("))'^"y,(s)(Z,(s), (iM,.(s)) (2.14) 
By (D.l) and (D.2), we have 

2j[^e2/o™)'^-y.(.)(F(X(.),yi(.),Zi(.))-F(X(s),y,2(,),z2(.)))ds 
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^ e'^o'^^^i-))^-Y,{s){F{Xis),Y,Hs),Zl{s)) - FiX{s),Y,\s) 
+ 2 r e'Io^^M-))^^Y^(s){F{Xis),Y,\s),Zl{s)) - F(X{s^^^^^ 



+d2^ /^e2/o'^(^("))'^"|Z^.(s)|2ds. (2.15) 
Choosing c = 26d2, we obtain from (j2.15p 

\e-fodmu))duY^^^^^2 ^ _ _i.)^^e-2/o'^W"))'^"(^(X(s))Z,(s),Z,(s))ds 

< g2/o^d(x(«))d«|y^^^^|2 _2 /"^d(X(s))e-2/o'^(^H)'^«y^(s)(Z^(s),dM^(s)) (2.16) 



Taking expectation on both sides of the above inequaUty and letting T tend to infinity, we 
obtain that 

We conclude that Y^{t) = Y^{t) and hence from ([2T6]) . Zl{t) = Z^{t). 
Existence: 

First of all, the assumption (|2.1L)p implies (see [?]) 

poo 

dLg] < oo. 

X Jo 

1°: There exists {pxi't),Qx{t)) such that 

dp,{t) = efo^^^(-))'^-^X{t))dLt+ < q,{t),dM.,{t) >, (2.17) 

and eio ^(^(«))^"p^(t) ^ as t ^ oo in L2(J7). 
In fact, let 

/oo 
e^o^^M-))d-<^{X{s))dL,\Ft] 

POO 

e^o 9(^H)rf«$(X(s))L, - / e^o 9(^H)'^«$(X(s))dL,|J-t]. 

(2.18) 

By the martingale representation theorem in [22], there exists a process qx{t)^ such that 

/"CxD /'OO 

ei'o '^(^("»'^"$(X(s))(iL,|J-i] = - e^o'ii^(^))'i^<^{X{s))dLs] 
Jo Jo 

+ / <q^{s),dM^{s)> . (2.19) 
JO 
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Then {px,qx) satisfies the equation (|2.17p . 
Moreover, 



eIo^W-))^-^X{s))dLs\Tt] 



Jo 

t f'°° 

= _e/o QiX{u))du-^^^ / q{X{u+t))du^j^(^^ ^ t))dLs+t\J^t] 

Jo 

t f'°° I - 

= -e-loiW^))'i-^Ex{t)[ e^o9(^H)<i«$(X(0)di.i] (2.20) 

Jo 

The last equality follows from the fact that Lt+g = Lt + Ls o Of. Therefore, 



t f °° t 

sup \px{t)\ < e^O <i{X{n))du 1^(^)1 . / g/„ g(X(«))dn^_^^]^ 

xgD xeD Jo 

Set M = sup |«>(x)| • sup / «'(^("))'^"(iLf]. 

In view of ([230]) . we have lim e/o('^+9)(^H)'='« = in L2(J1). 
Hence, 

g/o*d{x{n))d«^^^^) <^g/o*(d+g){x{«))d«^Q as t^oo, in L^{n). (2.21) 
2° : Set g{t,y,z) = F{X{t),px{t) +y,qx + z). Then 

(yi - y2){g{t,yi,z) - g{t,y2,z)) 

= {yi - y2){F{X{t),px{t) +yi,qx + z)- F{X{t),px{t) + ys, Qx + z)) 

< -di{Xmyi-y2\^. (2.22) 



and 



\git,y,zi)-g{t,y,Z2)\ = \F{X{t),pxit) + y,qx + zi) - F{Xit),pxit) + y,qx + zi)\ 

< d2\zi-Z2\. (2.23) 



Moreover, 



Ex[ / e^/o '^(^("))'^"|^7(X(t), 0, 0)|2dt] 
Jo 

< e2/o'^(^(«))'^«|F(X(t),p,(t),g,(t))|2dt] 

JO 

/•oo ^ 

< Ex[ e2/o'i(^W)'^"|E:(t)|2(;t] 

JO 

< oo. (2.24) 
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g satisfies all the conditions of the Lemma 2.2. Hence, there exist processes {kx,lx) such that 
dkx{t) = -g{t,kx{t)Mt))dt+ < Ut),dMx{t) >, 

and 

as t — )■ cxD. 

Putting Yx{t) = Pxit) + kx{t) and Zx{t) = qx{t) + Ixit), we find that {Yx{t), Zx{t)) satisfies the 
following equation 

dYxit) = efo^^^(^))'^^ct>{X{t))dLt - F{t,Yx{t),Zx{t))dt+ < Zx{t),dMx > . 

and 

Urn e/o'^(^("))'='"yi = 0. 



Corollary 2.1 Suppose all the assumptions in Theorem 2.1 hold. If, in addition, 

fOQ ^ 

supEx[ / eio '^(^("))'^"|i^(t)|2dt] < oo, 

x Jo 



then it follows that 



sup |y^.(0)| < oo. 

x€D 



Proof. 

As shown in the proof of Theorem 2.1, Yx{t) has the decomposition: Yx{t) = Px{t) + kx{t). 
Setting t = in ^a?M . it foUows that 



/■oo J 

|p.(0)| < Ex^t)[\ / e/o5-(^W)'^"ci>(x(/))dL;| 
Jo 

< ||<l>||oosupii;,[/ e^o^W-))'iudLi] 

X Jo 



< OO. (2.25) 

By Ito's formula, we obtain 

de^ /o d{X{u))du I |2 ^ _2g2 (^)^(^^ ^ 

+ 26^/0 ^(^("))^"A;^(t)d(X(t))dt + 2e2/o'^(^H)'^«A;^(t) < lx{t),dMx{t) > 
+ /o dixiu))du (^A{X{t))lx (t), /x.(t))di 

Choosing two positive numbers 6i and 82 such that > ^ and 5i+52 < S, similar calculations 
as in the proof of Theorem 2.1 yield that, for any t < T, 

£;,[e2/o'^(^(«))'^«|fe,(t)|2] + i(i _ -^)Ex[ Te^/o 

A 2Adi Jt 
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Setting t = 0, we have 



^020(2 -'0 



Let T — 00 to obtain that 

X \ZO2a2 X Jo J 

where the fact that e^^o '^(^("))'^"/c^(T) ^ as T — 7- cx), has been used. Hence, we have 
sup|ya;(0)| < sup|p^(0)| +sup|/c^.(0)| < 00. 

3 Linear PDEs 

Set 

L2 = - V • (^V) + 6 • V + g 

where h = {bi,...,bd) is a i?*^- valued Borel measurable function, and g is a Borel measurable 
function on R'^ such that: 

lD{\b\' + \q\)eL^{D), p>^. 
In this section, we solve the following linear boundary value problem: 

(3.1) 

where F and (b are bounded measurable functions on D. 



iy • (AVn)(x) + h ■ Vu{x) + q{x)u{x) = F{x), on D 
\^{x) = 4> ondD, 



It is well known that operator L2 defined on a bounded domain D with Neumann bound- 

du I 



ary condition = is associated with the quadratic form: 



£{f,9)- = - L2f{x)g{x)dx 



D 

{AVf,'Vg)dx— / b ■ V f{x)g{x)dx — / q{x)f{x)g{x)dx 

D Jd Jd 



Definition 3.1 A bounded continuous function u{x) defined on D is a weak solution of the 
problem ([3?!]) if u G W'^''^{D), and for any g £ C°°(D), 



£{u,g) = / (j){x)g{x)a{dx) — 1 F{x)g{x)dx, 
JdD Jd 

where a denotes the d — 1 dimensional Lebesgue measure on dD. 
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Consider the operator 



Lo = • (AVu) (3.2) 



on domain D with boundary condition = on dD. 

Lq is associated with a reflecting diffusion process {X^,P^). By [13], has the foflowing 
decomposition: 

dX^ = a(X?)dWt + ^VA{Xf)dt + 7iX?)dLl 



Lt - I I{xo€dD}dLg, (3.3) 







where the matrix (7(x) is the positive definite symmetric square root of the matrix A(x) and 
{^t}t>0 is a d-dimensional standard Brownian motion. 

It is well known that operator Lq is associated with the regular Dirichlet form: 



1 f dix du 



and the domain of 8^ is W^''^{D) := {u G L'^{D) : ^ G L^{D)}. 

The following lemma can be proved similarly as the Corollary 3.8 in [llj using the heat kernel 
estimates in 1211. 



Lemma 3.1 There exists a constant K > 0, such that 

sup^°[L°] < KVi and inf > 0. 

Moreover, we have sup^g^j E^[{L^)^] < K^t^ , for some constant > 0. 
Set MO = /(J a{X^)dWs and 

where b* is the transpose of the row vector b. 

The proof of the following two lemmas are inspired by that of the Lemma 2.1 and Theorem 
2.2 in dl]. 

Lemma 3.2 For t > 0, there are two strictly positive functions Mi{t) and M2{t) such that, 
for any x GD, Mi{t) < E^[Jq Z^dL^] < M2{t). Furthermore, M2{t) ^ as t ^ 0. 
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Proof. 
1°: Put 



M(t) = eIo<^-'HX%dM^>-y^';bA-^*{xO)ds ^ 
xGD jo 



(3.5) 



Then we have 



ZctdL' 



Ol 



x€D Jo 

< supii;°[max |M(s)|2]i .supii;°[e2|g|(t)(L0)2]5 

x€D 



XGD 0^^^* 



< supii;0[|M(t)|2]2.supii;0[e4u(t)]5-supii;°[(L°)^]4 (3.6) 



x&D 



xeD 



(///) 



By Khash'Minskh's lemma and Theorem 2.1 in [15], (/) and (//) are bounded if t belongs to 
a bounded interval. Because of S°[(L^)''] < K^t^ , we see that M2{t) := K{I){II)y/t is the 
required upper bound. 
2°: Since 

El[L^,f < El[ f M-\s)e-,{s)dL',] ■ E^,[ f M{s)e,{s)dL% (3.7) 
Jo Jo 



we obtain 



Here 



E'. 



M-\t) 



M{s)eq{s)dL°,] > 



E'xiL't? 



El[JlM-\s)e.,{s)dL^] 



(3.8) 



(3.9) 



By the proof of the first part, replacing Mt, q by A^^ and bA ^h* — q respectively, it is seen that 
there exists K{t) > such that sup^g;^ ^"[/^ M~^(s)e_q(s)(iL°] < K{t). 



As inf -Tj£'^[L^] > 0, we complete the proof of the lemma by setting Mi{t) = — 



Set G(x) := E^^if^^ Z,dL% 



Lemma 3.3 // there is a point xq G D, such that G{xo) < oo, then there are two positive 
constants K and /3 such that svlP^^J) E^[Zt] < Ke~^^ . 
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Proof. 

By Girsanov Theorem and Feymann-Kac formula, L2 = • (AV) + 6 • V + g is associated 
with the semigroup {Tt]t>Q, where Ttf{x) = El[Ztf{X^)] for / G L'^{D). 
By the upper and lower bound estimates of the heat kernel p2{t,x,y) associated with Tt in 
[21| . the following inequality holds, 

/ f{x)dx < < c / f{x)dx, (3.10) 

Jd J d 

where c is a positive constant. Since 

G{x) =Y,El[ZnEl,[ / ZsL\ds)]] > Mi(l) 

n=0 n=0 

and G(xo) < 00, there is a positive integer number A'' such that 
This implies 



j^El[ZN-i]m{dx)<^^. 



Thus 



sup^;°[Z7v] = sup^°[Zi£;^jZjv_i]] < c /" ^°[Z^_i]m((ix) < ^. (3.11) 



For any t > 0, there exists a positive number n such that £ [n — l,ra). Then by (j3.1ip . it 
follows that 

< :^^2[^t-^(n-i)] < ( sup El{zA^^ 

^ \x€D,0<t<N J ^ 



< 2 sup E^[Zt]e-'^Kn (3.12) 

x€D,0<t<N 



Theorem 3.1 If there exists xq € D such that G{xq) < 00, then there exists a unique bounded 
continuous weak solution of the problem (j3.ip ." 

Proof. 

Existence : 

Due to Theorem 3.2 in [1], there exists a unique, bounded, continuous weak solution U2 of the 
following problem: 



L2U2{x) = 0, on D 
l^ix) = <P ondD. 



(3.13) 
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Thus by the hnearity of the problem p.ip . we only need to show that the following problem 
has a bounded continuous weak solution: 



(3.14) 



L2Ui{x) = F{x), on D 
^ix)=0 ondD 

The semigroup associated with operator L2 is {Tt,t > 0}. By Lemma [331 we have 

sup |rtF(x)| = sup \E^,[ZtF{Xf)]\ < i^e-^*||F||oo. 

xeD x&D 

Then 

l-OO 

ui{x) := / TtF{x)dt 
Jo 

is well defined and has the following bound: 

sup \ui{x)\ < — ||F||oo. 
xgD P 

The function ui{x) is also continuous on D. 

In fact, fixing any x G D and e > 0, we can firstly choose a constant to > 0, such that 
sup2g£) I Jq" TsF{z)ds\ < |. And because TtgUi{x) is continuous, there exists a constant 6 > 0, 
such that for any y with |y — a;| < 6, \TtgUi{x) — TtQUi{y)\ < |. 
We find that 

/•oo 

TtUlix) = E',[ZtUiiX?)] = El[Zt E^o[ZsF{X',)]ds] 

Jo 

= [ E^^[Zt+sUiiX?^,)]ds 
Jo 



TsF{x)ds 

= ui{x)- I TsF{x)ds. (3.15) 
Jo 

For any y satisfying ly — x| < it follows that 

\ui{x) - ni(y)| < \Tt,m{x) - Tt,ui{y)\ + | /" ° TsF{x)ds\ + \ T T,F{y)ds\ < e. (3.16) 

Jo Jo 

This implies that the function ui is continuous on domain D. 

Denote the resolvents associated with operator L2 by {Gp,f3 > 0}. Note that 



POO 

G(sui{x) = / e-'^%ui{x)dt 
Jo 



00 />oo 

e-^*ui{x)dt- / e"^* / TsF{x)dsdt 
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We have 
Therefore, 



= -ui{x)- / / e'^^TsF{x)dsdt 
P Jo Jo 

-1 roc f'OC 

= -mix) - / TsF{x){ / e-^'dt)ds 

P JO Js 

= luiix)-^G^F{x). (3.17) 
I3{ui{x) - f3G^ui{x)) = l3GisF{x). 



hm / I3{ui{x) — PGi3Ui{x))ui{x)dx = Um / PG i3F{x)u{x)dx = I F{x)ui{x)dx < co. 
I^^o^Jd f^^°°JD Jd 

This imphes that ui G D{£) (see [16]) and ui is a weak solution of equation (j3.14p . By the 
linearity, u = ui + U2 is a bounded continuous weak solution of equation (jS.ip . 
Uniqueness : 

Let vi and V2 be two bounded continuous weak solutions of the equation (|3.1|) . Then vi — V2 
is the solution of equation ()3.13p with cp = 0. Then by the uniqueness of the equation ()3.13p 
proved in we know that vi = V2- □ 

4 Semilinear PDEs 

Recall that 



i,j=l j=l 



_d_ 

dxi 



and L2 = Li + q are two operators both defined on the domain D and equipped with the 
Neumann boundary condition ^ = on dD. 

(0, X(t), Pj, X £ D) is the reflecting diffusion process associated with the operator Li with 
the decomposition introduced in ()2.ip . 



In this section, we solve the following semilinear boundary value problem: 

L2u{x) = —G{x,u{x),Vu{x)), on D 



i|H(x) = ,/.(x) ondD ^^-^^ 



Let £{■, •) be the quadratic form associated with the operator L2: 

£{u,v) = — I < A\7u,\7v > dx — I < b,\7u > vdx — f quvdx. 
2 Jd Jd Jd 
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Definition 4.1 A bounded continuous function u{x) defined on D is called a weak solution of 
the equation (|4.ip if u € W^'^{D), and for any g G C°°{D), 

£{u,g)= / 4>{x)g{x)a{dx) + / G{x,u{x),Vu{x))g{x)dx. 
JdD Jd 

Recall that Lt is the boundary local time of X{t) defined in (j2.ip and is the boundary local 
time ofX^ in ([331) • 

As a consequence of the Girsanov theorem, we have: 

Lemma 4.1 Suppose that the function f satisfies Ex[Jq e-^o fi^('^))'^^(lL^^ < oo. Then it holds 
that 

Jo Jo 

where Mt was defined in (j3.5p . 

The following lemma is deduced from Theorem 3.2 in [3]. 

Lemma 4.2 Suppose that the function q G L^{D) and p > |. // there exists some point 
xq G D, such that 

Jo 

then it holds that 

sup^^[ / 9(^(«))'^«(iLt] < oo. 

X Jo 



Let G{x, y, z) : x R x R'^ — )■ i? be a bounded Borel measurable function. Introduce the 
following conditions: 

(H.l) (yi - y2){G{x,yi,z) - G{x,y2,z)) < -hi{x)\yi - y2p, 
(H.2) \G{x,y,zi) - G{x,y,Z2)\ < h2\zi - Z2\. 

Set h{t) = -hi{X{t)) + 5hl + q{X{t)) and h{t) = -hi{X{t)) + 5hl for some constant & > ^■ 

Theorem 4.1 Suppose that the conditions (H.l) and (H.2) are satisfied. Assume 

ExA e2/o(9(^H)+'^(«))«'«dt] < oo, for some xi£D, (4.3) 
Jo 

and there exists some point xq G D, such that 

Exoi / e^oi(^(^))d^dLt] < oo. (4.4) 
Jo 

Then the semilinear Neumann boundary value problem ()4.ip has a unique continuous weak 
solution. 
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Proof. 

Set 

Then 



and 



Note that 



G{X{t), y, z) := e^o e" /o e" /o lix {u))dt ^-^ _ 

(yi - y2)(G'(X(t), yi, z) - G{X{t),y2, z)) < -hi{x)\yi - (4.5) 
\G{X{t),y,zi) - G{X{t),y,Z2)\ < h2\zi - Z2\. (4.6) 



G{X{t),y,z) < e^oiW-))dt\\G\\^, 
By Theorem 2.1 there exists a unique process {Yx,Zx) satisfying 

dY,{t) = -GiX{t)X{t)J.it))dt + efo<im-))<i-<P{Xit))dm 
eIo'^(^)'i^Yx{t) ^0 as t ^ oo. 

Furthermore, Corollary 12.11 implies that sup 1^(0) < oo. 
From Ito's formula, it follows that 

d(e~/o9(^(«))'^*i;(t)) 
= -qiX{t))e-^o^'^^(^))dtYx{t)dt - e-^oi(^(^))dtG{X{t),Yx{t),Zx{t))dt 
+cP{X{t))dLt+ < e~fo'i(^(-))dtz,{t),dMxit) > . 

Setting Yx{t) := e" /o and Zx{t) := e" ''(x("))'^*Z^(t), we obtain 

dYx{t) = -{q{X{t))Yx{t) + G{X{t),Yx{t),Zx{t)))dt + <l){X{t))dLt+ < Zx{t),dMx{t) > . 

Moreover, 

g/o* h(u)dtY^^^^ = g/o* Ku)dt^- /,* q(x(u))dtY^^^^ = g/,f fe(x{«))dt-5> ^ Q t ^ oo. (4.7) 

So by Ito's formula, we have that, for any t <T, 

= eIohi-)duY^(^T)+ r efoh^-)d-(GiXxis),Yxis),Zx{t))+qiXxis))Yxis))ds 



/ e-^^^^^)'^^(t){X{s))dLs- h{s)e^o Hn)duY^^s)ds 
Jt Jt 



T 



e 



!^Ku)du^2x{t),dMx{t)). (4.8) 
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Put uo(x) = y^.(0) and voix) = Z^(0). 

Since 1^(0) = 1^(0), we know that uq is a bounded function on domain D. By the Markov 
property of X and the uniqueness of {Y^, Z^) , it is easy to see that 

Y,{t) = uo{X{t)), Z,{t) = vo{X{t)). 

So that sup ^ ll^ollco < co- 

xGD,t>0 

Now consider the following problem: 



L2u{x) = —G{x,uo{x),vo{x)), on D 
=</)(x) ondD 



(4.9) 



By Theorem 3.1, problem (|4.9p has a unique continuous weak solution u{x). Next we will show 
that u = uq. 

Since u belongs to the domain of the Dirichlet form associated with the process X(t), it follows 
from the Pukushima's decomposition that: 

du{X(t)) 

= -[G{X{t),uoiX{t)),vo{X{t))) + q{X{t))u{Xmdt + (t>{X{t))dL{t) + {Vu{X{t)),dMS)) 
= -[G{X{t),Y,{t),Z^{t))+q{X{t))u{Xm+(t>{X{t))dL{t) + {Vu{X{t)),dM,{t)) 

From the condition (14. 3p and the boundedness of u{x), it follows that 

lim E^[e2/o < ||n||^ lim E^[e^ ■foCh+<i)(n)du^ = q. 

t—^oo t—^oo 

By Ito's formula, it follows that, for any t <T, 
e/o'^(«)'^«u(X(t)) 

= elo M«)^«^(x(T)) + 1^ e/o '^H'^-[G(X(s), + qiXis))uiX{s))]ds 

T 

e/o '^H'^^ {Vu{X (t) ),dM^{t)). (4.10) 

Jt 
Set 

v^{t) = u{X{t)) - Y^{t) and R^{t) = Vu{X{t)) - Z^{t). 
Subtracting the equations ()4.8p from ()4.10p . we obtain the following equation: for any t <T, 

e/o^(")'^"-i;(X(t)) 

= e-fo '^(")'^"z;(X(r)) + / iq{Xiu)) - h{u))efo '^W'^"t;(X(s))(is 



oo 
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T 



Set g{t) = e-f^o ^^'^^'^'^v{t). Taking conditional expectation on both sides of (j4.1ip . we find that 
git) = EMT) - h{s)9{s)ds\Ft] 

= EMT){1- r Ks)ds) + r r h{s)h{si)g{si)dsids\Ft] 

Jt Jt Js 

= S.b(T)(l - his)ds + h{s)dsf ) 



l-T i-T i-T 

\3 



+ (-1)'^ / / / h{s)h{sx)h{s2)g{,S2)ds2dsxds\Tt\. 

Jt Js Jsi 



Keeping iterating, we obtain 



nl 



1 

k=0 

+ r r r :■ r h{s)h{si)..:h{sMSn)dSn...dSids\J't] 

Jt Js J Si J Sn~l 

Since E^[\g{T)\e^t \h\is)dsi^ < oo, letting n — )• oo, by dominated convergence theorem, it follows 
that 

Then 

v{t) = £;,[t;(r)e/t^('*(^)-'^W)'^^|J-t] < (holloo + ||n||oo)i^x[e^*''^(''(^))'l-^t]. (4.11) 
Hence, it follows that 

< e^oi(^(^))d^\v{t)\ < (ll^/olloo + llnlloo) lim E,[e^o i(^(^))d^\Tt]. (4.12) 
Since the condition ()4.4p implies 

lim EJe^o 'i(xis))ds^ ^ q 

T-)-oo 

we deduce that E^[efo 'i(^('^'))'^'^\v(t)\] = and hence v{t) = 0, - a.s.. 

Therefore, for any t > 0, we have u{X{t)) = Yx{t) and Vu{X{t)) = Zx{t) by the uniqueness 
of the Doob- Meyer decomposition of semimartingales. In particular, u{x) = Ex:[u{Xx{0))] = 
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-£'x[ia'(0)] = uo{x). This shows that u{x) is a weak solution of the equation (j4.ip . 

If u is another solution of the problem (j4.ip . Then the processes Yx{t) := u{X(t)) and Zx{t) := 

'Vu{X{t)) satisfy the following equation 

dY,{t) = -GiX{t),Y,it),Z,it))dt - 4>{X{t))dLt+ < Z,{t),dM,{t) > . (4.13) 

Set Y^{t) = ei'o«(^("))'^"y^(t) and Z^{t) = e/o 
By chain rule, it follows that 

dY^it) = -G{X{t),Y^{t),Z,{t))dt + efo'iW-))<i^^{X{t))dm + {Z^{t),dM,{t)) 
Moreover, because u is bounded, we have 

lim e^oh(^)d^Y^(t) = lim e^* ''(")'^"n(X(t)) = 0. 
Therefore, from the uniqueness of the solution of the BSDE in Theorem 2.1, we have 

Y.it) = Y^{t) ZS) = Z.it). 

In particular, 

u{x) = E^[Y^{t)] = E^[Y^{t)] = u{x). 



5 Semilinear Elliptic PDEs with Singular Coefficients 

Recall the operator 

L = iv • (AV) + B-V -V ■ {B-) + Q 

on the domain D equipped with the mixed boundary condition on dD: 

1 <9u , - , , , 
-—-{B,n)u{x) = 0. 

The quadratic form associated with L is given by: 

1 f du dv sr^ /" n / \ 



dv 



[ Bi{x)— — u{x)dx — [ Q(x)u(x)v{x)dx, 
,Jd oxi Jd 



where (•.•) stands for the inner product in L'^{D). 
The domain of the quadratic form is 

V{Q) = W^'\D) := {n : n G L\D), ^ G L\D),i= 1, ...,4. 
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Let {St, t > 0} denote the semigroup generated by L. 



In this section, our aim is to solve the fohowing equation: 



Lf{x) = -F{x,f{x)), onD 

<B,n> (x)/(x) = ^(x) on dD ^ ' 



Definition 5.1 A bounded continuous function f{x) defined on D is called a weak solution of 
the equation (fSTT]) if f ^ W'^''^, and for any g E C°°(_D), 

Q{u,g)= / ^{x)g{x)a{dx) + / F{x,u{x))g{x)dx. 

JdD Jd 

Here the function F : i?^ x — t- i? is a bounded measurable function and satisfies the following 
condition: 

(E.l) {yi-y2){F{x,y,)-F{x,y2)) < -ri(x)|yi - yaP- 



Recall the following regular Dirichlet form 

(5.2) 



associated with the operator Lq = \V{AV) equipped with the Neumann boundary condition 
1^ = on dD. 

The associated reflecting diffusion process is denoted by {fi, 5^^, X^^, 6*^, 7^, P^}. Here 6^ and 
7^ are the shift and reverse operators defined by 

X^M{^)) = Xf+.(a;),s,t>0 
X0(7?M) = XlM.s<t. 

The process {Xi)t>Q has the decomposition in (13 .Sp . The martingale part of is = 
jla{X^,)dWs. 

The following probabilistic representation of semigroup St was proved in [5] 

Stf{x) = E',[f{Xf) exp( [\A-'BnX^)dM^ + ( fiA'^ Bf {X^^)dM^^) o 7° 

JQ Jo 



i f\B-B)A-\B~Br{X^^)ds+ fQ{X^,)ds)] 
^ Jo Jo 



(5.3) 



E^ denotes the expectation under P^. 
Set 



Zt = exp(/ iA-^Bnx',)dM^ + {[ {A'^ B)* iX^,)dM^) o 
Jo Jo 

[\B-B)A-\B-Bnx^)ds+ fQ{X',)ds). 
^ Jo Jo 



(5.4) 



26 



By [3] and [21], there exists a bounded, continuous functions v G W^'^{D) satisfying that 
Jo 

Vv{X'^,)dMs + v{Xf) - v{X^) - / {A-^By{X'^^)dMs (5.5) 
Moreover, v satisfies the fohowing equations: for g E C^{D), 

/ < A\7v,\7g > {x)dx = / < B,\7g > {x)dx. (5.6) 
Thus the representation of St becomes: 

ft 

Stf{x) = e-^(^)^°[/(Xf)e^(^")exp( / {A-\B - B - AVv)ydM^ 

Jo 

-- [ {B - 13 - AVv)*A-\B - B - AVv)iX^)ds 

+ ^ (Q + ^(V7;)^(V^;)* -{B-B, Vv))iX'^)ds)] 
= e-^^^^St[fe^]{x). (5.7) 

Here, setting b := B - B - (AVv) and q := Q + ^{Vv)A(yv)* - {B - B, Vv), we see that St 
is the semigroup generated by the fohowing operator: 

L2 = ^V-{AV) + {B-B-{AVv))-V + {Q + ^(yv)A{Vvy-{B-B,Vv)) 

= • (^V) + b-V + q 

equipped with the boundary condition ^ = 0. 

In this section, we wih stick to this particular choice of b and q. 

Recall that 

M{t) = efo '^''^(^s)dMO-y^bA-^b*ixO)ds 
and set Zt = M{t)e^o i(^s)ds ^ 

Then from ([53]), it follows that Z{t) = Zte^(^?)-^(^o). 

Recall the operator Li = • (^V) + b ■ V with Neumann boundary condition, which is 
associated with the reflecting diffusion (X(t),P^). It is known from [T3] that 

dP^\^, = MtdP^lT,, 

and 

X{t) =x + [ a{X{s))dWs + [ il-VA + b){X{s))ds + [ -/{X{s))dL„ - a.s. 
Jo Jo 2 Jo 
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where {Wt} is a d-dimensional Brownian motion and Lt is the local time satisfying that Lt = 

jll9D{X{s))dLs. 

Lemma 5.1 Assume that there exists xq G D, such that 

/•oo 

Eli |Zi|2e/o(2Q-4n)(x{J)d«^^0] (5_8) 
Jo 

Then there exists a positive number e > 0, if \\B\\lp < e, the following inequality holds: 

supE^[ e2i'o(--i+9)(^("))<i«(ft] < oo. (5.9) 
xgd Jo 

Proof. 

= E°[Z'(i)e'^o(-2ri+9)(X(«))dMj 

< (7^£;0[z(t)e-2/o(''i(^("))'^"e^o('3+f<^v^-2(s--B),Vi,>)(xS)d«] 

< CiS°[Z2(t)e2/o(Q-2ri)(XS)dnji . ^0[g/o*<AVi,-2(B-J3),Vi,>(XO)d«]i 

By Lemma 3.3 and condition (15. 8p . there exists two constant C2,/3 > such that 

supii;,[Z2(t)e2/oW--i)(^S)<^«] < C2e-^*. 

Moreover, for p > by the Theorem 2.1 in [ihi, there exist two positive constants C3 and C4 

S IJ. ell t tl 3jt 

where C4 = c|| < AVv - 2{B - B),Vv > ||^p/2. 

Since IVwI^p < C|i?|^p(£)) (see [21j), there exists e > 0, such that l-Bl^^p^^)) < e implies C4 < /3. 
Thus dSSD holds. □ 

Theorem 5.1 Assume ()5.8p and for some point xq £ D 

/"OO 

ZsdL°]<oo (5.10) 

JO 

T/ien i/iere exisfo e > such that if \\B\\lp < e, the problem ()5.ip has a unique, bounded, 
continuous weak solution u{x). 

Proof. 

Existence: Set F{x,y) = e'"^^'' F{x,e~'"^^^y) and (f>{x) = e'"^^^^{x). 
From the boundedness of f , -F is also bounded. 
And F satisfies 

(yi - y2){F{x,yi) - F{x,y2)) < -ri{x)\yi - 2/2p- 
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Moreover, there is a constant c > 0, such that 

/"CXD /"CO 

oo>i^2o[/ Z^dL's] = Ell Z.e^(^^')-(x«)^^0] 

^0 ^0 

pQQ poo 

> cEli ZsdL^^]=cEl[ M,e/o''(^S)^«dLO] (5.11) 
Jo Jo 

By Lemma l4.lt we know that, at Xq G D, 

E.oi e/o^(^")'^«dL,] <oo. (5.12) 
Jo 

Furthermore, by Lemma |4.2| it follows that 

X Jo 

By Lemma |5.H the following condition is satisfied : 

Ecci e^^o(i-^^)(^(^))'i^dt]<oo, (5.14) 
Jo 

So F satisfies all of the conditions in Theorem 14.11 replacing G by F. Thus the following 
problem 



L2u{x) = —F{x,u{x)), on D 

Idu 
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has a unique bounded continuous weak solution u{x). 

Set f{x) = e~^^^^'u(x). Then we claim the function f{x) is the weak solution of the equation 

m- 

Because function v is continuous and bounded, f{x) is also continuous. From the fact that 
function u is the weak solution of the problem (j5.15p . we obtain, for any function ^ G C°°{D), 

£{u, e-^'V) = - / < ^Vn, V(e~''V) > - <b,Vu> e'"^ - e'^'qu^dx 

e-''(j)ipda+ / F(x,ti(x))e->(ix. (5.16) 
dD Jd 

As in the proof of Theorem 5.1 in [22], we can show that the left side of the equation (j5.16p 
equals to 

Q(/, ^) = l f [< AVf, V^j> - <B,Vu>^- <B,V^> f - QfiP]dx. 
^ Jd 

At the same time, by the definition of the function (\) and the right side of the equation 
()5.16p equals to 



/ ^i)da^ \ F{x,f{x))ipdx. 
JdD Jd 
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Thus it follows that, for any G C°°{D), 

Q{f,iP)= [ ^iPda+ I F{xJ{x))^dx. 

JdD J D 

which proves that function / is a weak solution of the problem (j5.ip . 
Uniqueness: 

If / is another solution of the problem (j5.ip . then u := e" f can be shown to be the solution 
of the equation (|5.15p . Then by the uniqueness of the problem (|5.15p proved in the Theorem 
14.11 we find u = u. Therefore, f = f.d 

6 solutions of the BSDE and Semilinear PDEs 

Recall the operator 




on the domian D equipped with the Neumann boundary condition ^ = 0, on dD. 

And {^l, Tt, X{t), Px, X G D) is the reflecting diffusion process associated with the generator 

Li. 

Then the process X(t) has the following decomposition: 

X{t) = X{0)+M{t)+ I b{X{s))ds+ [ An{X{s))dLs, P^-a.s.. 

Jo Jo 

Here b = {bi, ...,bd} with ftj = | J2j -^(^) -^t square integrable continuous 

martingale additive functional. 

In this section, we will consider the solutions of the BSDEs in Section 2 and use this 
result to solve the nonlinear elliptic partial differential equation with the mixed boundary con- 
dition. 

Let / : X i?+ X — 7- i? be progressively measurable. Consider the following conditions: 

(1.1) {y — y'){f{t,y) — f{t,y')) < d{t)\y — y'\'^ , where d(t) is a progressively measurable process; 

(1.2) E[/~e/o '^(«)'^"|/(s,0)|ds] < oo; 

(1.3) Px — a.s., for any t > 0, y — )■ f{t,y) is continuous; 

(1.4) Vr > 0, T > 0, Mt) ■■= sup \f{t,y)-f(t,0)\ G L\[0,T] x x dP^). 

\y\<r 

The following lemma is deduced from Corollary 2.3 in [2\. 

Lemma 6.1 Suppose a pair of progressively measurable processes (Y, Z) with values in Rx 
such that t ^ Zf belongs to L^([0,T]) and t — )• f(t,Yt) belongs to L^([0, T]), Px — a.s.. 
If 

Yt = (+ r fir, Y,)dr - T < Z^^dM^ >, (6.1) 
Jt Jt 
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then the following inequality holds, for < t < u < T , 

\Yt\<\Yu\+ Ysfis,Y,)ds- t{Zr,dMr). 
Jt Jt 

where y = ^I^y^oy 

The following lemma can be proved by modifying the proof of Proposition 6.4 in [2]. 

Lemma 6.2 Assume that conditions (I.l)-(1.4) with d{t) = 0. Then there exists a unique 
solution {Y, Z) of the BSDE 

Yt= [ f{r,Yr)dr- [ {Zr,dMr), for t<T. (6.2) 
Jt Jt 

Moreover, for each /3 £ (0, 1), Elsupt^j^ \Yt\'^] + E[{Jf^ |Z,.pdr) 2 ] < 00. 
Suppose /3 G (0, 1). 

denotes the set of real-valued, adapted and continuous process {lt}i>o such that 

\\Yf := E[sup\Ytf] < 00. 
t>o 

It is known that || • ||^ deduces a complete metric on . 

denotes the set of i?''- valued predictable processes {Zt} such that 



00 

2 - • 



z\\mp-=e[{ \Ztydt)2]<^. 







is also a complete metric space with the distance deduced by || • \\mP- 

Lemma 6.3 Under the same assumption as the Lemma 16. 2^ there exists a unique solution 
(y, Z) of the BSDE 

Yt = YT+ I fir,Yr)dr- [ (Z„dM,,), any t<T; 
Jt Jt 

limyt = 0, P-a.s.. (6.3) 

t— >oo 

Proof. Existence: 

By the Lemma 6.2 above, there exists {Y^, Z") such that, for < t < n, 

Y^= \ fir,Y,^)dr- / (Z,",dM,), 
Jt Jt 

and Yl" = = 0, for t > n. 
Fix t > and t < n < n + i, then 

rn+i j-n+i rn+i 

yn+i _ yn ^ / (^(^^ ^n+i) _ ^(^^ ^n))^^ _ / ^(^n+i _ ^n)^ ^^^^ ^ / ^(^^ g)^^ 

J t J t Jn 



31 



Set F"(r, y) = /{r, y + F^) - /(r, + /(r, 0)I{,>„}, = ^t"+' " and = - Zf. 

Then {yf,z'f) is the solution of the following BSDE: 

/•n+i j-n+i 

y^= F{r,y^)dr- {z^.dM,). (6.4) 



Jt Jt 
So that by the condition (I.l) with d{t) = 0, it follows from Lemma |6. II that 

rn+i rn+i 

\y7\ < / {y':,F^{r,y^))dr- {rr^z^Mr) 
Jt Jt 

rn+i J f'U+i 

< / ^f^(y^/(r,y- + y,")-/(r,y,"))dr+ / \fis,0)\ds 

Jt lUrl Jn 

n+i 

{y'^^zyMr) 

< / \f(s,0)\ds- {y^^z^Mr). (6.5) 

Jn Jt 

Taking conditional expectation on both side of the inequality, we got 

rn+i 



pn+r 

\y^\<E[ \f{sMds\Ft]:=M^, 

J n 



where M" is a martingale. Then by Doob's inequality and condition (1.2), it follows that, for 
/3g (0,1), 

1 f-n+i 

E[^nMP]<E[^nmtf\ < T^Ei \f{s,0)\dsf 

t t i- p Jn 

— )• 0, as n ^ oo. (6.6) 

Therefore, {Y^} is a Cauchy sequence under the norm || • ||^. So that there is a process Y 

such that E[supt \Yt - Yt^'f] -+ 0. 

This also implies that If — )• 0, as t — )• oo, Px — a.s.. 

Moreover, by the equation (16. 4|) . Ito's formula and the condition (I.l), it follows that 

/n+i 
{A{X{r))z^,z^)dT 

/•n+i j-n+i 

= 2 / (y-, F-(r, y^))dr - 2 / {y^ , z'^dM,) 
Jt Jt 

/•n+i rn+i 

< 2 {yl^,r{r,0))dr + 2\ {y^,z^dMr)\ 

Jn J t 

rn+i rOD 

< sup|y,"|2 + (/ |/(r,0)|dr)2 + 2| / (y,",<dM,)|, 

r Jn Jt 

and thus that 

rn+i g rn+i rn+i „ 

( / K\^dr)2 < ci[sup|y,"|'' + ( / \f{r,0)\dr)^ + | / < y,",z-dM, > 1 2]. 

Jt r Jn Jt 
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Taking expectation on both sides of the inequahty and applying the BDG inequahty, we obtain 



n+i fn+i , 

|y;?|2|z;^|2dr)4 



rn+t 

E[{ |z,"|2dr)f] 

/•n+i !• 

< c,{E[snp\y^\^]+E[{ \f{r,0)\drf])+C2E[{ 

■r Jt Jt 

< ciiE[snp\y:^f]+E[{ \f{r,0)\drf])+C2E[{sup\y';\^ Iz^^l'drf^] 

r Jn r Jt 

< (ci + |)(i?[sup|y,"n+i?[(y^ \f^r,0)\drf]) + -E[{J^ 



Therefore, we know that there is a constant C > 0, such that 



roo i-n+i 

E[{ \z^\'ds)2] < CE[snp\y^f + { \f{s,0)\ds)^] 

Jo t Jn 

pn+i 

< CE[snp\y^f] + CE[ \f{s,0)\ds] 

t Jn 



— )■ as n — )■ oo. 

So that {ZJ^} is a Cauchy sequence in M^. Let Z denote the hmit of {Z"}. 
At last, by the condition (1.3), we find that 

r f{t,Yndt^ r f{t,Yt)dt, P,-a.s.. (6.7) 
Jo Jo 

Therefore, (Y,Z) is the solution satisfies the BSDE (|6.3p . 
Uniqueness: 

Consider {Y,Z) and {Y',Z') are two solutions to (|6.3|) . Then by the same method as in the 
proof of Lemma 2.1, we can show that, 

Vt>0, \Yt-Yt\=0, P-a.s.. □ 

(1.5) The process d{t) is a progressively measurable process satisfying 

d{-) e L^[[0,T] X n,dt(g) P], for any T > 0. 



Theorem 6.1 Assume the conditions (I.l)-(1.4)- Then there exists a unique process {Y,Z) 
such that, 

Yt = YT+ [ f{r,Yr)dr- [ < Zr,dMr >, for any t<T; 
Jt Jt 

lim e^o d{u)duY^ = 0, P- a.s. (6.8) 

t— >oo 
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Proof. 
Existence: 

Set = e/o '^(")'^"/(t,e-i'o - d{t)y. Then 

(1) iy-y')if{t,y)-f{t,y'))<0; 

(2) /(t,0) =e/o'^(«)'^"/(t,0). Soi?[/;°|/(.,0)|ds] =S[/o°°e/o'^W'^"|/(t,0)|ds] <oo. 

(3) sup I / (t , y ) - / (t , 0) I < ^/^r (t) + 1 d(t ) I r , where the process i{jr{t)+\d{t)\r £ L^{[0,T]xn,dt^P), 
\y\<r 

for T > 0. 

Therefore, / satisfies ah the conditions of the Lemma 16.31 So there exists a pair of processes 
(Y, Z) satisfying the equation: 

Yt = YT+ I fir, Yr)dr- [ (Z, , dMr) , 
Jt Jt 

and obviously hm Yj = 0. 

By the chain rule and the definition of the function /, it follows that 

de- -fo dHduy^ ^ _ g- /J diu)duY^^^^ ^ ^g- d{u)du^^^ ^^^^ ^ 

Set Yt = e-/o'^(«)«'«yi and Zt = e" /o Then the process {Y,Z) is the solution to the 

equation ()6.8p . 

Uniqueness: 

The uniqueness of the solution to (|6.8p follows from the uniqueness of the solution to equation 
(lOD . □ 

Let G{x, y) : x R ^ R he a bounded Borel measurable function. Consider the following 
conditions: 

(H.l)' {yi -y2){G{x,yi,z) - G{x,y2,z)) < -hi{x)\yi -y2p, where hi G LP{D) for p > f . 
(H.2) y — )• G{x,y) is continuous. 

Theorem 6.2 Assume the Conditions {H.l)' and {H.2)' and that 
such that 

roo 

Jo 

Then the semilinear Neumann boundary value problem 

( L2u{x) = —G{x,u{x)), on D 
{ ^{x) = (t>{x) on dD 

has a unique continuous weak solution. 



there is some point xq G D, 
(6.9) 

(6.10) 
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Proof. 
Step 1 

Set G{X{t),y) = e/o <i{X{u))dtQ^^^ ^- q{X{u))dtyy rj,^^^ ^^^^^ ^^-^^g ^ unique solution Z^) 

to the following BSDE: 

for any T > and < t < T, 

%{t) = %{T)+ r G{X{t)X{s))ds- r e^o'ii^i-))dt^i^x{s))dLs 
Jt Jt 

and 

lim e' /o hi{Xiu))duY^ ^ Q ^^^^ 

t—>-oo 

The uniqueness follows from the uniqueness proved in Theorem 6.1. Only the existence of 
solution (Yx, Zx) needs to be proved: 

(a) Similarly as the proof of Theorem 2.1, we can show that there exists (Pxit), Qxit)) such that 

dpxit) = efoi(^i-))<i-^{X{t))dLt+ < qx{t),dMx{t) >, 

as t^oo, Px-a.s.. (6.11) 

(b) Set g{x,y) = G{x,y + Pxit))- Then it follows that 

{y - v'){9{x, y) - g{x, y')) < -hi{x)\y - y'f. 
The condition (|6.9p and Lemma 3.3 imply, for x G D, 

Ex[ / e^o(-'^^+'i)(^(^))<^^ds] < oo. 
Jo 

Furthermore, as the function G is bounded, we see that condition (/.2) is satisfied: 

roo 

Ex[ / e-/o '^i(^("»'^"|5(X(s),0)|ds] 

POO 

= Ex[ e'fo'^^(^(-'»d-\G(X{s),Pxis))\ds] 
Jo 

POO 

= Ex[ e-^oi-hi+i)ixiu))du^Q(^x{s),e--^oii^i^))'i^p^{s))\ds] 
Jo 

< \\G\\ooEx[ e^o(-'^^+i)(^i^))'i^dt] 

Jo 

< oo. (6.12) 

Obviously condition (1.3) is satisfied, i.e., y — g{x,y) is continuous. 
Moreover, the condition (1.4) is also satisfied. In fact, for any r > 0, 

Mt) =snp\G{X{t),y) -G{X{t),0)\ <2\\G\\^e^o'i(^t)dt^ 
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and for any T > 0, by the fact that q S U'{D) with p > ^ and Theorem 2.1 in |15j . 
< oo. 

Therefore, the function g{x,y) satisfies all of the conditions of Theorem 16.11 . There exists a 
pair of processes {yx{t), Zx(t)) such that for any T > and < t < T, 



Vxit) = vAT) + g{X{s),y,is))ds- {z^s) , dM^s)) (6.13) 

and 

lim e-^oh'(^('-))^^y^(t) = Px- a.s. (6.14) 



t^oo 



Put Yx{t) = Px{t) + yxit) and Zx{t) = qx{t) + Zx{t). It follows that {Yx{t), Zx{t)) satisfies the 
following equation 

dYxit) = e^o^'^^(-))^^cl){X{t))dLt - G{t,Yx{t))dt+ < Zx{t),dMx >, 
lim e- hiixiu))duY^ ^ Q 

t—>-oo 

Step 2. 

Put y^(t) := e-i'o9(^(«))'i*i;(t) and Z^(t) := e' Jo li^i^))'^^ Zx{t), we have 

dy,(t) = -F{X{t),Yx{t)) + cl){X{t))dLt+ < Zx{t),dMxit) >, 
where F{x,y) = q{x)y + G{x,y). Moreover, 

(,!o{-hi+q)(X{u))duY^{l^^ = g/o*(-hi+g)(X{M))(n)dtg-/o*g{X(n))dty^(^^^) 

Put ^0(2;) = ^^(0) and vq{x) = Zx{^)- 

Now as in the proof of Theorem, 4.1, we can solve the following equation 



00. 



L2u{x) = —G{x,uo{x)), on D 
1 ( 
2 



'if (a;) = Hx) on dD 



(6.15) 

and prove that the solution u coincides with uo{x). This completes the proof of the theorem. □ 



Suppose that F : R'^ x R R is a bounded measurable function and ri G U'[D). Con- 
sider the following conditions : 

(E.l) {yi- y2){G{x,yi,z) - G{x,y2,z)) < -ri(x)|yi - y2p; 
(E.3) y — )■ F{x, y) is continuous; 

Now, after establishing Theorem 6.2, following the same proof as that of Theorem 5.1, we 
finally have the following main result. 
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Theorem 6.3 Suppose that the function F satisfies the condition (E.l) and (E.2), and there 
exists xq E D such that 



poo 

<[/ Z,dLO]<oo. (6.16) 
Jo 

Then the following problem 

J Lu{x) = —F{x,u{x)), on D 



\ ^^{x)- < B,n> {x)u{x) = ^{x) ondD 
has a unique, bounded, continuous weak solution. 
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